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EARLY STAR FORMATION TRACED BY THE HIGHEST REDSHIFT QUASARS 
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ABSTRACT 

The iron abundance relative to a-elements in the circumnuclear region of quasars is regarded as a clock of the 
star formation history and, more specifically, of the enrichment by type la supernovae (SNIa). We investigate the 
iron abundance in a sample of 22 quasars in the redshift range 3.0<z<6.4 by measuring their rest frame UV Fell 
bump, which is shifted into the near-IR, and by comparing it with the Mgll A2798 flux. The observations were 
performed with a device that can obtain near-IR spectra in the range 0.8-2.4/im in one shot, thereby enabling an 
optimal removal of the continuum underlying the Fell bump. We detect iron in all quasars including the highest 
redshift (z=6.4) quasar currently known. The uniform observational technique and the wide redshift range allows 
a reliable study of the trend of the Fell/Mgll ratio with redshift. We find the Fell/Mgll ratio is nearly constant 
at all redshifts, although there is marginal evidence for a higher Fell/Mgll ratio in the quasars at z^6. If the 
Fell/Mgll ratio reflects the Fe/a abundance, this result suggests that the z^6 quasars have already undergone 
a major episode of iron enrichment. We discuss the possible implications of this finding for the star formation 
history at z>6. We also detect a population of weak iron emitters at z~4.5, which are possibly hosted in systems 
that evolved more slowly. Alternatively, the trend of the Fell/Mgll ratio at high redshift may reflect significantly 
different physical conditions of the circumnuclear gas in such high redshift quasars. 

Subject headings: quasars: emission lines — galaxies: evolution, high-redshift 



1. INTRODUCTION 

The iron abundance relative to a-elements has been regarded 
as a clock indicator of the star formation history (e.g. Hamann 
& Ferland 1999). Indeed most of the iron in local galaxies is 
thought to be produced by type la SNe, while a-elements are 
predominantly produced by type II SNe. The difference be- 
tween the progenitor masses of SNIa and SNII (hence in their 
life times) translates into an enrichment delay between the Fe 
and a-elements, which was generally thought to be about 1 Gyr 
(Wheeler, Sneden, & Truran 1989). However, more recent stud- 
ies have shown that SNIa are also enhanced in young stellar 
systems (e.g. Mannucci et al. 2003). In particular the SNIa 
progenitors may be as massive as 8 Mq (Greggio & Renzini 
1983), which have a lifetime of ^ 3 x lO^yr. As a consequence 
the maximum of iron enrichment may occur on relatively short 
timescales depending on the star formation history and on the 
initial mass function (Matteucci & Recchi 2001). 

The prominent emission lines observed in quasars allow us to 
investigate the metallicity of galactic nuclei at large cosmologi- 
cal distances (Hamann & Ferland 1999). The main iron feature 
in quasar spectra is the UV bump at 2200-3090A, which is the 
blend of several thousand Fell (and some Felll) lines. The ratio 
of the Fell (UV) bump to the Mgll A2798 doublet is sensitive 
to the Fe/a abundance ratio, but is also sensitive to other phys- 
ical properties of the emitting region such as microturbulence, 
density, and ionization parameter (Verner et al. 2003). As a con- 
sequence it is not simple to quantitatively derive the Fe/a ratio 
from the FeII(UV)/MgII ratio alone. The FeII(UV)/MgII ratio 
can be used as a relative indicator of Fe/a among quasars, un- 
der the assumption that the physical properties of the emitting 
regions are the same. 

At z>3 the UV iron bump is shifted in the infrared. Sev- 



eral authors have measured the Fell/Mgll ratio in high redshift 
quasars by means of near-IR spectroscopy (Aoki, Murayama, & 
Denda 2002; Dietrich et al. 2002, 2003; FreudHng, Corbin, & 
Korista 2003; Iwamuro et al. 2002; Thompson, Hill, & Elston 
1999). Yet, the main difficulty in several of these studies (with 
the exception of Dietrich et al. 2002, 2003) is that most infrared 
instruments can only cover a fraction of the near-IR spectrum in 
a single observation, and the broad Fell pseudo continuum can 
well cover the entire observed spectral range, making it very 
difficult to recover the true underlying continuum. 

In this letter the results of the analysis of near-IR spectra of 
22 high redshift quasars will be presented. These spectra were 
obtained with an instrumental setup that can cover the full near- 
IR spectrum from 0.8/im to 2.4/im in one shot (and with very 
high throughput), therefore allowing a good estimate of the con- 
tinuum underlying the Fell bump. Our observations provide a 
relatively large sample of Fell/Mgll measurements over a wide 
redshift range and obtained with the same technique and fitting 
method. Therefore these data allow for a reliable investigation 
of the redshift dependence of the Fell/Mgll ratio. Through- 
out the paper we adopt Hq = 70 km s"'Mpc"', J^m = 0.3 and 
f7A = 0.7. 
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Fig. 1. — Near infrared spectrum of the highest redshift quasar 
SDSSJl 14816.64+525150.3 at z=6.4. The upper thick line shows the observed 
spectrum, where regions of bad atmospheric transmission have been omitted, 
while the overimposed thin solid line shows the best fit to the continuum. Fell 
and Mgll, as discussed in the text. The dotted line shows the power-law com- 
ponent (¥x ocA",a = — 1.9), while the dashed fine shows the total continuum 
fit (i.e. by including the Balmer continuum). The lower part of the figure shows 
the single components of the fit and, more specifically, Fell bump, Mgll line 
and Balmer continuum. 

2. OBSERVATIONS 

The observations were obtained with the Near Infrared Cam- 
era Spectrograph (NICS) at the ItaHan Telescopio Nazionale 
Galileo, a 3.56 m telescope. Among the various imaging and 
spectroscopic observing modes (Baffa et al. 2001), NICS of- 
fers a unique, high sensitivity, low resolution observing mode, 
which uses an Amici prism as a dispersing element (Oliva 
2003). In this mode it is possible to obtain the spectrum from 
0.8/um to 2.4/im in one shot. The throughputof the Amici prism 
is nearly two times higher than other more commonly used dis- 
persers. The spectral resolution with a 0".75 slit, as it was in 
our case, is 75 (i.e. 4000 km s"') and nearly constant over the 
whole wavelength range. This observing mode is appropriate 
to study the near-infrared continuum of faint sources, and to 
detect broad (~ 5000 km s"') emission lines in faint quasars as 
well as pseudo-continua such as the Fell bump which require a 
careful subtraction of the underlying continuum. Observations 
were performed in three observing runs: 2002 November 7-9; 
2003 February 25-29; and 2003 May 23-26. Typical integration 
times ranged from ^20 minutes, for bright objects, to ^2 hours 
for faint objects. We include also the data of SDSSJl 044-0 125 
obtained previously at the same instrument and already pub- 
lished in MaioUno et al. (2001). Several quasars were observed 
more than once on different nights to check for any instrumental 
or observational artifacts in the individual spectra. Wavelength 
calibration was performed by using an argon lamp and the deep 
telluric absorption features. The telluric absorption was then 
removed by dividing the quasar spectrum by a reference star 
spectrum (FV-GV) observed at similar airmass. The intrinsic 
features and slope of the reference star were then removed by 
multiplying the corrected quasar spectrum by a spectrum of the 
same stellar type from the library by Pickles (1998), smoothed 
to our resolution. Absolute flux calibration was obtained by 
using the photometry on the acquisition image or through pho- 
tometry reported in the literature. 

^ The best redshift fitting the observed wavelength of MgH in our spectrum is 6. 
Barth et al. (2003). 



In total 22 quasars were observed. These were selected 
among the brightest (observable during each run) in each of the 
following redshift ranges: 3.0-3.6, 4.3-5.1, 5.8-6.4, which en- 
sure that the Mgll and a good fraction of the Fell are observed 
outside the deepest atmospheric absorption regions (with the 
only exception of SDSS1044-0125 at z=5.78 whose Mgll is in 
a bad atmospheric band, but which could be used to derive the 
average quasars spectrum, see next section). 
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Fig. 2. — Average quasar spectra (residuals after continuum subtraction) 
around the Fell UV bump obtained by grouping the objects in the three main 
redshift bins and by normalizing each spectrum to the continuum at the wave- 
length of Mgll. Beneath each panel we also show the number of individual 
quasar spectra which where used to create the average spectrum. Regions 
where the number of quasars is low correspond to the spectral segments more 
seriously affected by tiie strong atmospheric absorption in that redshift range 
(see text). 

3. DATA ANALYSIS AND FITTING TECHNIQUE 

In Fig.l we show the resulting spectrum of the highest red- 
shift quasar SDSSJl 148 16.64H-525 150.3, at z=6.40 (Fan et al. 
2003)^, along with the identification of the major spectral fea- 
tures. Some blue-shift is observed for the CIV line. This shift 
may be due to an imperfect correction of the atmospheric ab- 
sorption dip which occurs at the same wavelength; however, 
Barth et al. (2003) also find the same effect, suggesting that 
the shift may be real. The spectrum of the highest redshift 
quasar shows a prominent Fell bump. All the other spectra 
and the detailed analysis of their spectral features will be pre- 
sented in a forthcoming paper. Here we focus on the analy- 
sis of the UV iron complex and of the Mgll A2798A blend. 
The continuum fitting we performed was similar to the method 
adopted by Dietrich et al. (2002, 2003) and Wills, Netzer, & 

I, in agreement with the redshift obtained by WiUott, McLure, & Jarvis (2003) and 
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Wills (1985), i.e. we included both a powerlaw and a Balmer 
continuum. The continuum was fitted in the spectral regions 
free of emission features, and away from regions of very bad 
atmospheric transmission. In Fig.l we show the power-law 
(Fa cx A", a = -1.9) plus Balmer-continuum fit to the spec- 
trum of SDSSJ114816.64H-525150.3. Note that in the highest 
redshift quasars (z>5.8) there is a flattening in the continuum 
blueward of CIV, that is also observed in the optical spectra 
of quasars at 3<z<5 (Francis et al. 1991; Vanden Berk et al. 
2001), possibly as a consequence of intervening metal absorp- 
tion hnes. At variance with other authors, who prefer to fit 
the powerlaw by sampling the continuum redward of 3100A 
and the continuum between CIV and Lya (but leaving positive 
continuum residuals between CIV and Fell), we prefer to fit 
separately the continuum on either side of CIV with two dif- 
ferent power-laws. This also ensures that the continuum fit of 
the quasars at z>5.8 is consistent with the continuum fit of the 
lower redshift quasars. 

In Fig. 2 we show the averages of the quasars residual spec- 
tra after continuum subtraction (and normalized to the contin- 
uum flux near the Mgll line) and grouped in the three main red- 
shift ranges discussed above. In each spectrum the regions of 
strong atmospheric absorption were discarded, but the spread in 
redshift allows to cover continuously the whole spectral region 
around Fell andMgll. Beneath each average spectrum we show 
the number of quasar spectra which populate each spectral re- 
gion. While both the low redshift and the high redshift average 
spectra clearly show a prominent Fell bump, the average of the 
intermediate redshift quasars is characterized by a significantly 
lower intensity of the Fell bump (middle panel). This might 
be partly due to the region at A < 2500A being poorly popu- 
lated since severe atmospheric absorption separating the J and 
H bands affects several quasars in this spectral region at these 
redshifts. However, the Fell emission close to Mgll is not af- 
fected by this problem. Therefore we believe that the weakness 
of Fell in the average spectrum of intermediate redshift quasars 
is real. 

The Mgll line and the residuals in the region around the Fell 
bump (2200-3090A rest frame) were fitted, for each object, by 
using the Vestergaard & Wilkes (2001) iron template, smoothed 
to the velocity dispersion of Mgll, and a broad gaussian at the 
wavelength of the Mgll fine. These fitting components for the 
case of SDSSJ114816.64H-525150.3 are shown in Fig.l. The 
Fell flux was integrated between the rest-frame wavelengths 
2200 and 3090A. In Table I we report flie Fell/Mgll ratio for aU 
quasars observed by us. We also list the rest-frame luminosity 
XLx at 1450A (note that thanks to our wide spectral coverage 
we can measure this quantity directly in several quasars, or by 
performing only a minor extrapolation). 

We note that the continuum and Fe fitting techniques, as well 
as the instrumental setup, have often been different among pre- 
vious works on this subject. As a consequence, a comparison 
of the results between these different studies may be subject to 
strong uncertainties (although some authors attempt to correct 
as much as possible for these systemic differences, e.g. Dietrich 
et al. 2003). Within this context we emphasize that our survey 
provides uniform and homogeneous Fell/Mgll measurements 
for a set of 22 quasars spanning a wide redshift range (3.0-6.4) 
obtained with the same instrumental setup (which covers the 
entire near-IR in one spectrum), same spectral resolution (con- 
stant over the whole spectral range) and same fitting technique. 
Therefore they represent a self-consistent sample which is well 



suited to study the trend of Fell/Mgll with redshift. 

Some of the quasars observed by us were also observed 
by other authors and deserve some comparison, even with 
flie caveats discussed above. BR 2237-0607 and BR 0019- 
1522 were observed by Dietrich et al. (2003) who find higher 
Fell/Mgll ratios than us, but the inconsistency is not large 
when errors are taken into account. BR 2237-0607 was also 
observed by Iwamuro et al. (2002), who obtain a ratio con- 
sistent with ours. SDSSJ103027.10H-052455.0 was observed 
with HST/NICMOS by Freudhng, Corbin, & Korista (2003) 
and they derived a much lower value for the Fell/Mgll ratio 
(2. 1 ± 1 . 1) flian found by us (8.64 ± 2.47). However, Freudling 
et al. did not measure Mgll directly in this object, but instead 
estimated the Mgll flux from the measured CIII] flux and an as- 
sumed CIII]/MgII ratio similar to the average of other quasars. 
Also, they did not sample the continuum redward of the Fell 
bump. SDSSJ114816.64H-525150.3 was observed by Barth et 
al. (2003), who found a lower Fell/Mgll value (4.7 ±0.4) than 
ours, though marginally consistent with errors (~ 1. 5a). Their 
spectral coverage does not extend enough to sample the contin- 
uum outside the iron bump, this might have led them to overes- 
timate the underlying continuum. 

In Fig. 3 we show the average of the Fell/Mgll ratios within 
the three redshift ranges discussed above. In the same figure 
we also plot the Fell/Mgll ratio obtained by Dietrich et al. 
(2003) from composite optical spectra of quasars at lower red- 
shift matching nearly the same luminosity range as the objects 
in our sample. No trend of the Fell/Mgll with luminosity was 
found within the luminosity range of our sample, which spans 
about 1.5 orders of magnitude. 



Fig. 3. — Filled symbols show the average of Fell/Mgll ratios of the quasars 
in our three main redshift bins. Hollow symbols show the Fell/Mgll ratio ob- 
tained by Dietrich et al. (2003) from optical composite spectra of quasars at 
lower redshift. 

4. DISCUSSION 

Some evidence for iron super-solar abundance was found by 
Wills, Netzer, & Wills (1985) in low to intermediate redshift 
quasars. In our spectra we detect iron emission also in all high 
redshift quasars. The ratio of Fell/Mgll is nearly constant over 
the redshift range zero to 6.4. The analysis of the high redshift 
quasars indicate a possible increase in the Fell/Mgll ratio at 
redshifts z ~ 6. If the Fell/Mgll flux ratio is tracing the Fe/a 
abundance ratio, this result would imply that the iron abundance 
is nearly constant or even slightly increasing from redshift zero 
to 6.4, i.e. when the age of the universe was about 800 Myr. 
While this finding is difficult to reconcile with the scenarios ex- 
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pecting a delay of 1 Gyr or more for the iron enrichment, mod- 
els of rapid and strong star formation at high redshift offer a 
plausible explanation. Indeed, Matteucci & Recchi (2001) have 
shown that models of early formation of ellipticals, with very 
efficient star formation and lasting '^0.4 Gyr, imply a peak of 
SNIa production (hence of iron enrichment) already at 0.3 Gyr 
after the onset of star formation. They have also shown that in 
the case of an instantaneous burst the SNIa production peaks at 
only 50 Myr after the burst. As a consequence, our results can 
be explained with two possible scenarios: 

1. A major episode of star formation in these highest-z 
quasar hosts occurred at z > 9, to allow the estimated 
0.3 Gyr required for SNIa to enrich the ISM within the 
context of models for star formation in elliptical hosts. 
The assumption that high redshift quasars are hosted in 
young elliptical galaxies may be reasonable, based on 
the studies of quasar hosts at lower redshift (Kukula et 
al. 2001; Dunlop et al. 2003). Nonetheless, it is also 
important to keep in mind that there is no direct obser- 
vational evidence that these highest-z quasars reside in 
"normally" evolving elliptical galaxies. 

2. An "instantaneous" burst of star formation occurred as 
late as z ^ 7, and SNIa had enough time to enrich the 
ISM within 50 Myr of the burst. The ideal case of 
a single instantaneous burst is not likely to represent 
the scenario of spheroids and galaxies formation in gen- 
eral, for which a star formation extended in time is more 
plausible. However, a single instantaneous burst may 
have occurred in the nuclear region of these quasars, 
causing a rapid iron enrichment hmited mostly to the 
central region. 

It is interesting to note that the first scenario (early formation 
of ellipticals) implies a redshift for the first major episode of 



star formation (z>9) which is consistent with the redshift of the 
re-ionization (at z« 20 ± 10) inferred from the recent Wilkin- 
son Microwave Anisotropy Probe (WMAP) results (Bennett et 
al. 2003). While the second scenario would be more consis- 
tent with the secondary re-ionization (at Z'~6) proposed by Cen 
(2003) and Wyithe & Loeb (2003). 

We also find a population of quasars at z~4.3-5.1 with low 
iron emission (Table 1). These objects are responsible for the 
apparently lower Fell bump in the average spectrum of the 
quasars in this redshift bin (Fig.2). Other quasars with such 
low Fe emission in this redshift range were also found by Diet- 
rich et al. (2003) and by Iwamuro et al. (2002) (see also Fiore et 
al. 2003). Statistically we cannot make any strong statement on 
the fraction of quasars with low Fell. However, they could rep- 
resent a population of quasars hosted in stellar systems which 
evolved more slowly than the quasars with higher Fe emission. 

An alternative explanation for the Fell/Mgll trends discussed 
above is that the physical conditions of the Broad Line Re- 
gion in high redshift quasars are significantly different than in 
quasars at lower redshift. In this case the observed Fell/Mgll 
ratio at high redshift would reflect the combined effects of 
abundance variations and variations of other physical param- 
eters such as density, ionization parameter, and microturbu- 
lence. Disentangling these effects requires our observations to 
be complemented with the detection of the optical iron bump 
(Vemer et al. 2003), which is shifted beyond the K-band at z>5. 
Observations sensitive enough to detect the optical iron bump 
at these redshifts will probably have to await the James Webb 
Space Telescope (JWST). 

We thank the referee M. Dietrich for helpful comments. Part 
of this work was developed during the Guillermo Haro Work- 
shop in 2003. RM acknowledges partial financial support from 
the Italian Space Agency (ASI). YJ and RM acknowledge sup- 
port from CONACyT grant J32178-E. 
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Table 1 
FeII/MgII measurements 



Name z Fell/Mgll log ALa(1450A) 



BR 0019-1522 


4, 


.52 


2.02±1, 


.01 


46, 


.76 


PSS J0134+3307 


4, 


.53 


1.68±0 


.42 


46, 


.87 


SDSS J021 102.72-000910.3 


4, 


.73 


4.10±2 


.05 


46 


.17 


SDSS J075618.14+410408.6 


5, 


.09 


3.43±1 


.37 


46 


.54 


SDSS J103027.10H-052455.0 


6 


.28 


8.65±2, 


.47 


46 


.68 


SDSS J104433.04-012502.2 


5, 


.78 


a 




46. 


.88 


SDSS J104845.05H-463718.3 


6, 


.23 


8.03±1, 


.22 


46. 


.81 


SDSS Jl 14816.64-1-525150.3 


6, 


.40 


8.20±2 


.10 


46. 


.98 


SDSS J130608.26H-035626.3 


5, 


.99 


9.03±2 


.26 


47. 


.32 


GB 6 15007H-4848 


3, 


.20 


4.17±0, 


.63 


47. 


.63 


BRJ1603H-0721 


4, 


.38 


4.06±0 


.81 


46. 


.86 


SDSS J160501. 21-01 1220.6 


4, 


.92 


5.30±1 


.35 


46. 


.46 


PSS J1633H-1411 


4, 


.35 


3.23±1, 


.28 


47. 


.02 


HS 1649H-3905 


3, 


.05 


7.40±0, 


.74 


46. 


.72 


SDSS J173352.23H-540030.5 


3, 


.42 


7.74±1, 


.16 


46. 


.99 


SDSS J173744.88H-582829.6 


4, 


.94 


3.55±1 


.42 


46. 


.75 


S5 1759H-75 


3, 


.05 


5.01±2 


.01 


47. 


.41 


PSS J2154H-0335 


4, 


.36 


5.83±1, 


.75 


47. 


.81 


SDSS J220008.66H-001744.8 


4, 


.77 


5.50±1, 


.65 


46. 


.56 


BR 2237-0607 


4, 


.56 


3.00±1 


.20 


47. 


.42 


LBQS 2231-0015 


3, 


.02 


5.77±0 


.86 


47. 


.06 


SDSS J234625.67-001600.5 


3, 


.50 


11.7±2 


.92 


46. 


.90 



^The Mgll of this object was not measured because it is located in the 
atmospheric absorption between H and K. However this quasar was used 
to obtain the average spectrum in Figure 2. 



